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We review the recent theoretical developments towards understanding the Mott phases and quan- 
tum phase transition of extended Bose-Hubbard models on lattices in two spatial dimensions . We 
focuss on the description of these systems using the dual vortex picture and point out the crucial 
role played by the geometry of underlying lattices in determining the nature of both the Mott phases 
and the quantum phase transitions. We also briefly compare the results of dual vortex theory with 
quantum Monte Carlo results. 
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I. INTRODUCTION 



The superfluid to Mott insulators transitions of 
strongly correlated lattice bosons systems, described by 
extended Hubbard models, in two spatial dimensions 
have recently received a great deal of theoretical inter- 
est. One of the reasons for this renewed attention is 
the possibility of experimental realization of such mod- 
els using cold atoms trapped in optical lattices^ 2 -. How- 
ever, such transitions are also of interest from a purely 
theoretical point of view, since they provide us with a 
test bed for exploring the recently developed theoretical 
paradigm of non-Landau-Ginzburg- Wilson (LGW) phase 
transitions^. In the particular context of lattice bosons 
in two spatial dimensions, a general framework for such 
non-LGW transitions has been developed and applied to 
the case of square lattice^. The application of this frame- 
work to triangular—, kagome^ lattices has also been car- 
ried out. These theoretical works has also supplemented 
by quantum Monte Carlo (QMC) studies which serves 
as a rigorous test for predicted theoretical results. Such 
studies, for extended Bose Hubbard models, has been 
carried out for square^, triangulares, kagom o 11 ' 12 , and 
hexagonal lattices 13 . In this review, we present a brief 
overview of the theoretical ideas developed in the above- 
mentioned works focussing on role of the geometry of the 
underlying lattices in determining the nature of both the 
Mott phases and the quantum phase transitions. We also 
briefly compare the results of dual vortex theory with 
QMC results. 

The typical paradigm of non-LGW transitions that has 
been proposed in the context of lattice bosons in Refs. 
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Ulalaul is the following. It is proposed that for non-integer 
rational fillings / = p/q of bosons per unit cell of the un- 
derlying lattice (p and q are integers), the theory of phase 
transition from the superfluid to the Mott insulator state 
is described in terms of the vortices which are non-local 
topological excitations of the superfluid phase, living on 
the dual lattice. 14 : 15 These vortices are not the order pa- 
rameters of either superfluid or Mott insulating phases 
in the usual LGW sense. Thus the theory of the above 
mentioned phase transitions are not described in terms 



of the order parameters on cither side of the transition 
which is in contrast with the usual LGW paradigm of 
phase transitions. Also, as first explicitly demonstrated 
m Refs.H although these vortices are excitations of a fea- 
tureless superfluid phase, they exhibit a quantum order 
which depends on the filling fraction /. It is shown that 
the vortex fields describing the transition form multiplcts 
transforming under projective symmetry group (projec- 
tive representations of the space group of the underlying 
lattice and that this property of the vortices naturally 
and necessarily predicts broken translational symmetry 
of the Mott phase, where the vortices condense. Since 
this translational symmetry breaking is dependent on the 
symmetry group of the underlying lattice, geometry of 
the lattice naturally plays a key role in determining the 
competing ordered states of the Mott phase and in the 
theory of quantum phase transition between the Mott 
and the superfluid phases. 

In what follows, we shall develop the above-mentioned 
ideas concentrating on square and Kagome lattices. The 
relevant results for the triangular lattice can be found 
in Ref. [f| To this end, we consider the extended Bose- 
Hubbard Hamiltonian 
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Here t is the boson hopping amplitude between nearest 
neighbor sites, U is the on-site interaction, V denotes the 
strength of the nearest neighbor interaction between the 
bosons and ^ is the chemical potential. We note at this 
point that a simple Holstein-Primakoff transformations 
maps this boson model (Eq. [1]) , in the limit of hardcore 
bosons (U — > oo), to XXZ models with ferromagnetic 
J x and antifcrromagnetic J z interaction in a longitudinal 
magnetic field Bi 

h xxz = -j x J2(sfs* + sis]) + J^sfs; 

(ij) (ij) 
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where J x > and J z > are the strengths of transverse 
and longitudinal nearest neighbor interactions and B\ is a 
longitudinal magnetic field. The parameters J x , J z and 
Bi of Hxxz can be mapped to those of iJb oson • 'Jx — 
2t, J z = V, and Bi — zV (/ — |), where z denotes the 
coordination number of the underlying lattice and / is 
the average boson filling. 



II. COMPETING MOTT STATES AND 
QUANTUM PHASE TRANSITION 

The basic premise underlying the theories developed in 
Refs. lUIIlfilQ is that the quantum phase transition of the 
extended Bose-Hubbard model occurs due to destabiliza- 
tion of the superfluid phases by proliferation of vortices 
which are topological excitations of the superfluid phase. 
To analyze such a transition, one therefore needs to ob- 
tain an effective action for the vortex excitation. Such an 
action can be obtained by performing a duality analysis 
of the Bose-Hubbard model-' 1 '-. The final form of the 
dual action Sd can be written as^ 
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where ipt, are the vortex field living on the site b of the 
dual lattice, A^ is the U(l) dual gauge field so that 
e Tl/ \A„Ab\ = rii where n, is the physical boson den- 
sity at site i, J2 P denotes sum over elementary plaque- 
tte of the dual lattice, A M denotes lattice derivative along 
fx = x,y, t, and / is the average boson density. Note that 
the vortex action Sd is not self-dual to the boson action 
obtained from H^ oson . Therefore we can not, in general, 
obtain a mapping between the parameters of the two ac- 
tions, except for identifying e TV \l\ v A^\ as the physical 
boson density 4 . We shall therefore classify the phases of 
this action based on symmetry consideration and within 
the saddle point approximation as done in Refs. 

The transition from a superfluid ((ipb) = 0) to a 
Mott insulating phase in Sd can be obtained by tun- 
ing the parameter r. For r > 0, we are in the super- 
fluid phase. Note that the saddle point of the gauge 
fields Abp in action corresponds to t TV \lA v A\\ = /, so 
that the magnetic field seen by the vortices is pinned 
to the average boson filling /. Now as we approach 
the phase transition point r = 0, the fluctuations about 
this saddle point {(ipb) = 0, A\ y — fx) increase and 
ultimately destabilize the superfluid phase in favor of 
a Mott phase with (ipf,) 7^ 0. Clearly, in the above 
scenario, the most important fluctuations of the vortex 
field ipb are the ones which has the lowest energy. This 
prompts us to detect the minima of the vortex spec- 



trum by analyzing the kinetic energy term of the vor- 
tices i/kinctic = -VvJ2b,a (i't +a e 2 ' TlAb ipb + h.c) , where 
the sum over a is carried out over the sites of the dual 
lattice which are nearest neighbors to b. The analysis 
of -ffkinctic therefore amounts to solving the Hofstadtcr 
problem on the dual lattice which has been carried out in 
Refs. 13, [jj andlrj fl9ll2"ol for square, triangular and Kagome 
lattices respectively. 



A. Square lattice 

The dual lattice in this case is also a square lattice with 
its sites being at the center of the plaquettes of the di- 
rect lattice. Thus here one needs to solve the Hofstadter 
problem on a square lattice with / = p/q flux quanta 
per plaquette of the dual square lattice. It has been 
shown in Ref. that for / = p/q, there are q minima 
within the magnetic Brillouin zone at the wave-vectors 
Q = (k x a,k y a) — (0,2irlp/q), where / runs from to 
q — 1, k XjV are the wavevectors of the magnetic Brillouin 
zone and a is the lattice spacing. The eigenfunctions 
^(k x ,k y ) corresponding to these minima can be easily 
constructed^. 

Once the positions of these minima are located, the 
next task is to identify the low lying excitations around 
these minima which are going to play a leading role in 
destabilizing the superfluid phase. These excitations, 
represented by the bosonic fields <pi , has definite transfor- 
mation properties under the symmetry operations of the 
underlying dual lattice^. For the case of square lattice, 
the relevant symmetry operations are translations along 
x and y directions by one lattice spacing (T x and T y ), re- 
flections about x and y axes (I x and I y ) and rotation by 
7r/2 (R-n/2)- R has been shown in Ref. U that the trans- 
formation properties of cpi under these transformations 
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Finally, we need to construct an effective LGW action 
in terms of the low-energy fields (pi which is invariant 
under all the symmetry operations. For the sake of def- 
initeness, we shall focuss on the case q = 2 and q = 3 
from now on. First, let us consider the case q — 2 for 
which, the effective action reads 
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a 8 th order term Ls in the effective action since this is 
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FIG. 1: Plot of Sp(r) on the sites, links and plaquettes of the 
direct lattice: the lines represent links of the direct lattice. 
As discussed in the text, these values represent the boson 
density, kinetic energy, and ring-exchange amplitudes respec- 
tively. The choice between the states B and C is made by a 
eighth order term in the action. 

the lowest order term, invariant under all the symmetry 
operations, which determines the relative phase of the 
vortex fields go and Q\ . 

The Mott phases, which are stabilized for s < 0, 
can be obtained for q = 2 by a straightforward analy- 
sis of L4 and L$. To concisely present the results for 
these phases, it is convenient to define a generalized den- 
sity 5p(r) = Yjrn.n=-\.\Pmn^ [mr:cJrnrv) \ where p mn 

S (|Q mn |) uj mn / 2 J2i=o <P*i <Pi+nOJ lm denotes the most gen- 
eral gauge-invariant bilinear combinations of the vortex 
fields if 1 with appropriate transformation properties of 
square lattice space grourA As shown in Ref. 4 the val- 
ues of 6p(r) on sites of the dual lattice (integer r x and r y ) 
can be considered a measure of the ring-exchange ampli- 
tude of bosons around the plaquette whereas those with 
r half-odd-integer co-ordinates represent sites of the di- 
rect lattice, and the values of Sp(r) on such sites measure 
the boson density on these sites. Finally, r values with 
r x integer and r y half-odd-integer correspond to vertical 
links of the square lattice (and vice versa for horizontal 
links), and the values of Sp(r) on the links is a measure 
of the mean boson kinetic energy; if the bosons represent 
a spin system, this is a measure of the spin exchange 
energy. We also note two fundamental points from the 
definition of p m n- First, as long as (<fi) ^ for at least 
one non-zero Q mn , the Mott state is characterized by a 
non-trivial density wave order. Second, the relative phase 
of the boson fields gi which plays a crucial role in deter- 
mining the nature of the density wave order is fixed by 
the 8 th order term Lg. 

The Mott phases for q — 2 is shown in Fig. [T] We note 
that there are three possible phases. For 701 > 0, only 
one of the two vortex fields condenses and the resulting 
Mott state (state A in Fig. [T]) has a charge density wave 
order. For 701 < 0, both the vortex fields condense and 
this leads to two possible Mott states (states B and C 
in Fig. [T]). For these states, all sites of the direct lattice 



are equivalent. These therefore represent valence bond 
solid (VBS) states with columnar or plaquette VBS order 
parameters. 

Next, we address the issue of quantum phase transi- 
tion between the superfluid and the Mott state for q = 2. 
Such a transition, within the premise of dual vortex the- 
ory discussed here, is described by S v and not, as would 
be expected by LGW paradigm, by an action written in 
terms of order parameters fields at either side of the tran- 
sition. The phase transition can either be first or second 
order. If it turns out to be second order and if v < 0, 
the relative phase of the vortex fields, which is pinned 
in the Mott phase by the dangerously irrelevant (in the 
renormalization group sense) 8 th order term Lg, becomes 
a gapless mode at the transition 4 . Thus such a transi- 
tion has an emergent gapless mode at the critical point. 
It provides an example of a deconfined quantum criti- 
cal point^ and can be shown to accompanied by boson 
fractionalization 4 -. 

The Mott phases for q = 3 can also be similarly 
obtained 4 -. Here the quartic term in the effective La- 
grangian which determines the Mott states are given by 

U = ^(|^„| 2 + M 2 + |^| 2 ) 2 

TOl / * * 9 * * 2 * * * 2 

+ — VPoVi </>2 + PiVzPa + ViVof-PoPi + c - c - 
- 2|^lV| 2 -%ilV 2 | 2 -2|^| 2 |< A) | 2 ). (5) 

Note that here the ordering of the Mott states are com- 
pletely determined by £4 and retaining higher order 
terms are not essential. It turns out that there are two 
sixfold degenerate VBS states with columnar or diagonal 
VBS orders as shown in Fig. [2] The superfluid-Mott in- 
sulator quantum phase transition is again described by 
a dual vortex action, but does not lead to deconfined 
quantum criticality as pointed out in Ref. 0. 

QMC studies on Bose-Hubbard models on square lat- 
tice with nearest neighbor interaction and near half-filling 
has been carried out in Ref. 8. These studies indicate a 
checkerboard Mott state and a strong first order transi- 
tion between the superfluid and checkerboard Mott state. 
Analogous QMC studies on triangular lattice has also 
been carried out in Refs. 



B. Kagome lattice 

The dual lattice corresponding to Kagome is the well- 
known dice lattice which has three inequivalent sites de- 
noted commonly by A, B and G^. To obtain the min- 
ima of the vortex spectrum, we therefore need to solve 
the Hofstadter problem on the dice lattice. Here we shall 
concentrate only for / = 1/2 and 2/3 (or equivalently 
/ = l/3). 

For / = 1/2, it has been show n 19 i 20 that the vortex 
spectrum on dice lattice do not have well-defined min- 
ima and collapse to three degenerate bands. Physically, 



FIG. 2: Plot of 5p(r) for q = 3. 




FIG. 3: Aharanov-Bohm cages for vortices in a dice lattice. 
See text for details. 



the collapse of the vortex spectrum can be tied to local- 
ization of the vortex within the so-called Aharanov-Bohm 
cages as explicitly demonstrated in Ref. 113 . An example 
of such a cage is shown in Fig. [3] A vortex whose initial 
wavepacket is localized at the central (white) A site can 
never propagate beyond the black sites which form the 
border of the cage. This can be understood in terms of 
destructive Aharanov-Bohm interference: The vortex has 
two paths — > 1 — > 3 and — > 2 — > 3 to reach the rim 
site 3 from the starting site 0. The amplitudes from these 
paths destructively interfere for / = 1/2 to cancel each 
other. Thus the vortex remain within the cage. Such dy- 
namic localization of the vortex wavepackets, termed as 
Aharanov-Bohm caging, makes it energetically unfavor- 
able to condense vortices. Thus superfluidity persists for 
arbitrarily small values of t/V . In the language of spins, 
this also explains the absence of S x ordering for XXZ 
model in a Kagome lattice for Bi = and J z 3> J x ■ Such 
a persistence of superfluidity has also been confirmed by 
QMC studies^. 

In contrast to / = 1/2, the vortex spectrum 
has two well-defined minima for / = 2/3 within 
the magnetic Brillouin zone located at [k x a, k) = 
(0, 7r/3) and (2tt/3, 2tt/3), where k = V3k y a/2. Thus the 
low energy properties of the vortex system can be char- 



FIG. 4: Density wave state with 9 by 9 order (full pattern 
not shown) for / = 2/3. Up/down arrows denotes singly 
occupied/empty sites. See text for details. 

acterized in terms of the two low-energy bosonic fields <pi 
and (fi2 similar to the case of square lattice at / = 1/2. 
However, as shown in Ref. @, in contrast to the case of 
the square lattice, the space group of symmetry transfor- 
mations of the dice lattice involves translations by lat- 
tice vectors u = () and v = () (T u and T v ), rotation 
by 7r/3 (Rtt/3), and reflections around x and y (I x and 
Iy). The transformation properties of the vortex fields 
under these symmetry operations are the following: T u : 
ipi -> ipi exp(-i7r/3), ip 2 -> ^ 2 exp(i7r/3), T v : ipi — > 
V?iexp(i7r/3), cp 2 -> ip 2 exp(-wr/3), I x : yJi(2) ~> ¥>i(2)> 

h ■ Vl(2) ^2(1) ) and : ^1(2) -» ^(lJT- 

The simplest Landau-Ginzburg theory for the vortex 
fields which respects all the symmetries is£ 





= £(2) + L (4) + L (6) 
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Here we find that L v turns out to be the lowest-order 
term which breaks the U(l) symmetry associated with 
the relative phase of the vortex fields. 

A simple power counting shows that is marginal 
at tree level. Unfortunately, the relevance/irrelevance of 
such a term beyond the tree level, is not easily determined 
analytically^. If it so turns out that is irrelevant, the 
situation here will be identical to that of bosons on square 
lattice at / = 1/2. The relative phase of the vortices 
would emerge as a gapless low energy mode at the critical 
point. The quantum critical point would be deconfincd 
and shall be accompanied by boson fractionalization 4 . 
On the other hand, if is relevant, the relative phase 
degree of the bosons will always remain gapped and there 
will be no deconfinement at the quantum critical point. 



5 



Here there are two possibilities depending on the sign 
of u and v. If u,v < and w > 0, the transition may 
become weakly first order whereas for u > 0, it remains 
second order (but without any deconfincmcnt). QMC 
studies seem to indicate a weak first order transition in 
this case^i. However, it is possible that a second order 
quantum phase transition can still be possible at special 
points on the phase diagram and this issue remains to be 
settled! 1 - 1 ^ 3 .. 

The vortices condense for r < signifying the onset 
of Mott phases. It turns out that the Mott phase for 
v > do not lead to any density wave order—. For v < 0, 
there are two possible Mott phases as demonstrated in 
Ref. H One of these phases with 9 by 9 ordering pattern 
is shown in Fig. [4j Here the bosons are localized in red, 
blue, green (closed circle) and black (open circle) sites 
whereas the green (open circle) and black (closed circle) 
sites are vacant leading to a filling of / = 2/3. The net 
occupation of the hexagons takes values 0, 4 and 2 as 
shown in Fig. [U Note that interchanging the occupa- 
tions of all the black and green sites of this state (while 
leaving the red and the blue sites filled) has the effect 
of 2 <-> 4 for the boson occupation of the hexagons la- 
beled 2 and 4 in Fig. [4] while leaving those for hexagons 
labeled unchanged. The resultant state is degenerate 
(within mean-field theory) to the state shown in Fig. 2] 
This allows us to conjecture that the effect of quantum 
fluctuations on the state is shown in Fig. [4] might lead to 
stabilization of a partial resonating state with a 3 by 3 
R — 3 — 3 ordering pattern^. Such a state is indeed found 
to be the true ground state in exact diagonalization^ and 
QMC studies 11 . Note that the mean-field analysis of vor- 



tex theory which neglects quantum fluctuations can not 
directly give this partially resonating state. 



III. DISCUSSION 

In conclusion, we have briefly reviewed theoretical de- 
velopments on duality analysis of Bose-Hubbard and 
equivalently XXZ models on two dimensional lattices. 
We have focussed mainly at predictions such an anal- 
ysis for square and Kagome lattices at half and one-third 
filling. We have compared the predictions of such dual 
vortex theories with QMC results wherever possible. Fi- 
nally, we would like to point out a couple of issues that 
still remain unsettled in this regard. It is yet unknown, 
whether it is possible to have a second order quantum 
phase transition in these models as predicted by the vor- 
tex theory Direct QMC studies&SiiSiiiiiSiSS as well as 
several other studies 2 ^ seem to find at best a weak first 
order transition. However, it has been recently conjec- 
tured that at least for the Kagome lattice, a second order 
transition might occur at special particle-hole symmetric 
points in the phase diagrams 23 . Since QMC studies, at 
least for large enough system sizes, can only reach close 
to this point, this issue is not settled yet. Second, it re- 
mains to be seen whether such duality analysis for boson 
theories with U(l) symmetry can be extended to either 
models with Fermions and/or bosonic models with higher 
symmetry group. 
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S. Isakov, Y.B. Kim, R. Melko, S. Sachdev, and S. Wessel 
for collaborations on related earlier works, and K. Damlc, 
T. Senthil, and A. Vishwanath for helpful discussions. 
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